Quantum simulators, continuous-time automata, and translationally invariant systems. 



K. G. H. Vollbrecht 1 and J. I. Cirac 1 
1 Max-Planck Institut fur Quantenoptik, Hans-Kopfermann-Str. 1, Garching, D-85748, Germany 

(Dated: February 1, 2008) 

The general problem of finding the ground state energy of lattice Hamiltonians is known to be 
very hard, even for a quantum computer. We show here that this is the case even for translationally 
invariant systems. We also show that a quantum computer can be built in a ID chain with a fixed, 
translationally invariant Hamitonian consisting of nearest-neighbor interactions only. The result of 
the computation is obtained after a prescribed time with high probability. 

PACS numbers: 



The difficulty of simulating the dynamics of quantum 
systems by classical means was recognized by Feynman 
[l| more than two decades ago. He proposed to use an- 
other system, a quantum simulator, to overcome this 
problem, introducing several visionary ideas about quan- 
tum computation. At the same time, very powerful clas- 
sical techniques have been discovered, which allows us to 
tackle important problems in many-body physics. One 
of the most important questions in this context is to de- 
cide which kind of problems can (or cannot) be efficiently 
simulated by quantum (or classical) simulators. 

While it is clear that such a device can simulate the 
dynamics very efficiently Q, it seems that it cannot be 
used to prepare ground states of arbitrary nearest neigh- 
bor interacting Hamiltonians. In fact, very recently it 
has been shown that if that was possible even in one spa- 
tial dimension then a quantum computer would be 
able to solve all NP (and even QMA) problems efficiently, 
something which seems to be unreasonable. One may ar- 
gue that systems in Nature are not so general since they 
typically posses certain symmetries (like homogeneity, or, 
equivalently, translational invariance) which restrict very 
much the Hamiltonians we are interested in and thus it is 
not so surprising that some of them cannot be efficiently 
simulated. 

Another related question is whether a particular quan- 
tum simulator (eg. one that is translationally invariant 
(TI) and with nearest neighbor interactions only) may 
be as powerful as a general one [1] . This question can be 
answered in the positive if one shows that it can perform 
every quantum computation efficiently (with a polyno- 
mial overhead in terms of the number of qubits). In 
fact, we have previously shown that this is the case if 
one is able to change the evolution Hamiltonian with 
time 0, [||]. Very recently, Werner et al. 0| have also 
shown that it is possible by alternating two kind of dis- 
crete gates, a result which generalizes previous ones on 
cellular automata. However, it still remains to be seen 
if with a fixed Hamiltonian that is TI and which only 
includes nearest neighbor interactions it is still possible 
to perform arbitrary quantum computations efficiently. 

In this paper we address the two problems mentioned 
above. First, we show that by demanding that a ID 



Hamiltonian is TI, a quantum computer is not more ef- 
ficient in preparing its ground state. Note that in order 
to be able to accommodate the number of parameters 
which define QMA problems, we must relax the condi- 
tion of only nearest-neighbor interactions, although we 
still keep only three-body (sites) interactions. Thus, our 
results imply that the homogenity that typically appears 
in Nature is not enough to make it simulatable. Second, 
we show that it is possible to build a quantum computer 
based on a static, TI Hamiltonian with only nearest- 
neighbor interactions, as long as one is able to prepare 
arbitrary product states and measure the sites indepen- 
dently. This result extends previous ones on cellular au- 
tomata, but as opposed to those, has no simple classical 
analogue since we deal here with continuous time evolu- 
tion. The TI quantum computation scheme can be seen 
as a combination of a cellular automata approach with a 
continuous-time quantum walk [l|, 0] ■ 

Ground state energy in TI systems: We start out by 
showing that finding the ground state energy E of any 
TI Hamiltonian in a ID chain is very difficult (as we 
increase the number of sites), even for a quantum com- 
puter or simulator. In fact, we show that for each QMA- 
complete problem one can find a TI Hamiltonian, H n , 
such that: Eq = if the answer to the problem is 
'yes'; Eq > l/(Poly(n)), if the answer is 'no'. Here 
Poly(n) denotes any polynomial in n. This implies that 
any effective routine for finding the ground state energy 
would give us a possibility to efficiently solve all prob- 
lems in the complexity class QMA. Our result heavily 
relies on the recent discovery |3[ that QMA-complete 
problems can be encoded in the ground state energy of a 
Hamiltonian h n describing a ID chain of n system with 
d = 12 levels each and with nearest-neighbor interac- 
tions, such that its minimum eigenvalue, A m i n , is either 
zero or A m j n > Poly(n). Our strategy is to build H n out 
of h n by only increasing the dimension d by a factor of 
two. 

We take at every site one additional qubit and define 
the Hamiltonian H l n = \l)(l\i®h£'"" i+n) . That is, if the 
extra qubit at site i is in the state |1), then the Hamil- 
tonian h n is applied to the (i, . . . , i + n) particles in the 
chain in this order, where we identify the sites separated 
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by n + 1 sites. Furthermore, we define the Hamiltonian 
H' that only acts one the extra qubits as 

n(n — 1 

where k, k' , k" run from site f to n. Now, we take 

1 " 

H n = H' + - } H l n . 
n 

i—l 

This Hamiltonian is translationally invariant, and con- 
tains up to three-body interactions only. It is easy to 
show that it fulfills the conditions mentioned above, since 
its eigenvectors have the form \ki,k2, ■ ■ .) ® \4>), where 
|fci, k2, ■ ■ •} is a product state on the extra qubits only 
and \cf>) is a state on the rest of the chain. In case h n has a 
zero eigenvalue, with corresponding eigenvector \4>o), we 
can construct an eigenvector of H n , |1, 0, Q, . . . ,0) ® |0o), 
also with E n — 0. Otherwise, it is easy to show 
that the minimal eigenvalue of H n is lower bounded by 
min ( A " ti " , — tt^ — tt I , given that if there is no or more 
than one extra qubit in 1, this is penalized in energy by 
H', whereas if there is only one, then this is penalized by 

Continuous-time automata: Now, we introduce a pro- 
grammable quantum computation scheme for an infinite 
chain of quantum systems of dimension d — 30. The 
program is encoded into the initial state, while the time 
evolution is fixed and given by a universal TI Hamil- 
tonian with nearest-neighbor interaction only. We will 
start by introducing a simple standard quantum com- 
puting scheme. Later on, we will show how this quantum 
computer can be simulated by a continuous time evolu- 
tion. 

We consider a simple quantum computer with an n— 
qubit 'hard disk' and a read/ write head, that we call the 
pointer. This pointer can be moved such that we can 
address single qubits. Furthermore, the pointer has an 
internal quantum state, a qubit. To perform any quan- 
tum computation we write a program consisting of five 
different commands: L) The pointer is moved one site to 
the left; R) the pointer is moved one site to the right; 
S) the qubit at the position of the pointer and the inter- 
nal state of the pointer are swaped; G) a G-Gate Q on 
those two qubits is applied. The G-gate allows for arbi- 
trary quantum computations if it can be applied between 
any two qubits. In this simple model this can be accom- 
plished by loading qubits with the S-command into the 
internal pointer state, which can then be moved to any 
other qubit. We now encode this quantum computer into 
a higher dimensional chain in the following way: Every 
site in the chain has three registers. A 'qubit' register 
having dimension 2. It acts like the normal qubit of the 
quantum computer. A 'pointer' register having dimen- 
sion 3. Here we encode the pointer, where indicates no 
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FIG. 1: A 7-qubit quantum computer encoded into a one 
dimensional chain. The pointer is represented by the 1, the 
program is encoded right of the quantum computer. 

pointer and 1/2 the presence of a pointer with an internal 
qubit state. Finally, a 'program' register having dimen- 
sion 5: One for 'e=empty' and the rest for commands 
{L, R, S, G}. The total dimension of one site defined this 
way matches 30. 

We choose n neighboring sites to be the 'quantum com- 
puter', i.e., the qubit-registers of those sites correspond 
one to one to the qubits of the quantum-computer we 
want to simulate (see Fig. 1). The pointer registers are 
in state |0) everywhere, except for one site which contains 
|1). The program is written in the program register, an 
area to the right of the 'quantum computer', where the 
commands are arranged in the order they should be exe- 
cuted from left to right. The rest of the program registers 
are filled with |e). 

We assume now a TI Hamiltonian with nearest neigh- 
bor interactions, H = J^. Hi, where Hi denotes a two-site 
Hamiltonian acting only on the sites i and i + 1 and is 
given by 

H= ]T \e) i {C\®\C)H. 1 {e\®U? i+1 + h.c (1) 

Ce{L,R,S,G} 

The first two operators only act on the two program reg- 
isters at sites i, and i + 1, respectively, while U c, s are 
unitary operators acting both on the qubit and pointer 
register at those two sites. The unitaries U s and U G 
act only on the first site i. A swap resp. G-gate is ap- 
plied to both the i-th qubit and the internal state of the 
pointer in case a pointer is present, and the identity op- 
eration otherwise. U swaps the two pointer sites while 
U L swaps the two qubit sites. Let us assume now, that 
we apply one of those unitaries on every pair of sites 
(i,i + 1) starting from the right end of the chain going 
step by step to the left. II s and U G will do nothing un- 
til we reach the position where a pointer is present in 
the first of the two participating sites. In this case, the 
corresponding gates are applied. The U R gate will swap 
the two pointer states at every position. As a result the 
pointer, marked by |1), is moved one step to the right. 
Finally, U L will swap the two qubit registers. As a result 
the whole qubit-register chain is moved one step to the 
right. Relative to the 'quantum computer' position, the 
pointer is moved one step to the left. Therefore, if such 
a sequence of unitaries is applied from right to left the 
corresponding set of commands is applied to the 'quan- 
tum computer'. Note, that it is sufficient to start this 
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sequence of unitaries at any position to the right of the 
'quantum commuter' and to stop it once they are to its 
left. Our claim is now that this is exactly what the time 
evolution basically does. 

To this aim, we consider the time evolution operator 
e -tHt _ . Applied to the initial state \start) 

we just end up in a linear combination of states of the 
form H n \start). Since H is defined as a sum of terms 
of the form hf :— \e i} Cj + i)(C,, e i+ i\ <g> Uf i+1 (and the 
hcrmitian conjugates), we can write the result of a time 
evolution as a linear combination of states of the form 
h n . . . /13/12 hi\ start), where every h k e {hf,hf}. The 
effect of any such hf on a state is the following: it ei- 
ther moves the command C from place (i + 1) to i and 
applies U c to both sites + 1), or it maps the state 
to zero if no C is at i + 1 and no empty space |e) at 
i. In the same manner hf^ results in a state where the 
command C is moved from i to i + 1 while U c ^ is ap- 
plied. Therefore h n . . . hzh2h\\start) is either zero, or a 
state where several of the commands have moved while 
applying the assigned unitaries (or their conjugates) all 
the way from their initial to their final position. Note 
that such a command can only move, if the correspond- 
ing neighboring site is in the state |e). This implies that, 
up to some |e)'s in between, the order of the commands 
in the program register can never change. Further note 
that if a command moves to the left and afterwards to 
the right, we end up in the same state since U c and 
U c ^ cancel each other. This implies that two states 
h n . . . h 3 h 2 hi\start} , g m ■ ■ ■ 9392gi\start) are equal, iff the 
configuration of commands in the program registers are 
equal. 

Assume now that we measure the program registers 
in the standard product basis after some time evolution. 
According to the above discussion we can conclude just 
from the configuration of the program register what has 
happened to the pointer and qubit registers. In par- 
ticular, suppose that we measure a configuration where 
all the commands that were initially to the right of the 
'quantum computer'(see Fig. 1), are found to the left of 
the 'quantum computer'. In this case the whole program 
has been executed (in the right order) to the 'quantum 
computer' and we can read out the result. If this is not 
the case, we can in principle continue time evolution and 
repeat until we found a positive result. To boost the 
probability of success, we can increase the program code 
by some irrelevant code, e.g. by adding \L)'s that do not 
effect the result of the computation if carried out after 
the real program. This irrelevant code will act as a kind 
of barrier, that prevents the real program code to move 
to the right and forces it to do the computation. 

In order to verify that this is an effective way of car- 
rying out quantum computations, we have to guarantee 
that the probability of finding the whole program to the 
left of the 'quantum computer' is finite after a finite time 



(that only grows polynomially in the number of qubits). 
To this end we have to solve the above problem and calcu- 
late (a lower bound to) the probability of success. Note, 
that given the initial state, the whole system stays in a 
subspace that can be labelled by only the configuration 
of the commands in the program register. So, we can 
map our system onto a chain of qubits, where |0) means 
'empty' and |1) means 'command'. Note that, since the 
order of the commands stays unchanged under the evolu- 
tion, we do not have to distinguish between the different 
commands. The Hamiltonian simplifies to H — ^2 t Hi 
with Hi — |lj0j + i)(0jlj + i| + h.c. At this point we see 
that the efficiency of the computation does not depend 
on the specific program we want to run. Let us take 
this system as electrons in a lattice, where |1) stands 
for an 'electron' and |0) for an empty site and Hi is a 
hopping term, allowing the electron to hop from one site 
to a neighboring site. As usual for fcrmions, there are 
no hopping terms allowed resulting in two electrons sit- 
ting at the same site. Thus, we end up with a system 
of non-interacting fermions in second quantization. To 
solve it, it is simple to go back into first quantization, 
where we just have to consider a single electron Hamil- 
tonian, H = J2k \k + l)(k\+ h.c, and all the effects com- 
ing from the Pauli principle are automatically implied 
by the anti-symmetrization. The single electron problem 
can be easily solved. Let us assume M sites with peri- 
odic boundary conditions. The eigenvectors are given by 
ijj q = — = J2 X e l Or xq \x) with corresponding eigenvalues 

e(q) = 2co S (^q). 

Let us now look at M sites where the first N sites 
are filled with each one electron. Remember, that these 
N electrons correspond to the program and enough ir- 
relevant code to force the program to go into the de- 
sired direction. An electron sitting at site y is written 
in terms of the eigenvectors as <j> y = J2 q ' l Pq e ~ t ^' V9 = 
jj J2 q k e l v( x -y)i\ x ). After the time evolution, the state 

is changed into 4> y , t = jjY, q ,k et ^ (x ^ v)q+le(q)t \ x ) ■ Wc 
assume now N electrons sitting in the first N sites 
what leads to the state tpo — S[\<j>i, <f>2, (f>3, ■ ■ ■ , <Pn)], 
where S denotes the fermionic anti-symmetrization op- 
erator. After waiting a time t we end up in the state 
ipt = S[\<f>i tt ,<h,t,<l>3,t,---,<l>N,t)]- We now want to com- 
pute the probability p\ to find particle I at time t still 
in one of the first N sites. Note that, since the vectors 
<j) v j are still orthogonal, the reduction to the first parti- 
cle gives us an averaging over all the starting positions 
y. This probability yields 

N N 

J* = a?EE K<M*>I 2 - ( 2 ) 

y— 1 x—1 

In the Appendix it is shown that by choosing t to be pro- 
portional to N, e.g. t — 500(W, we can bound this prob- 
ability to be smaller than a fixed constant, e.g. p\ < 0.3. 
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So we can guarantee to find particle 1 after a polyno- 
mial time with probability bigger than p = 0.7 outside 
the starting area. From this we can calculate the ex- 
pected number of electrons that will be found in those 
areas. They will be Npi and iV(l — pi) = Np. For a suc- 
cessful computation a fixed number k of electrons have 
to leave their starting area. But what is the probability 
to find more than k electrons outside the starting area 
if we know that the expectation value is Np? Let us 
assume the worst case scenario: We either get (k — 1) 
electrons with probability pf or N electrons with proba- 
bility p s . Since we know that pj(k — 1) + p s N = Np we 
can conclude that the success probability p s > -j^r^r 
which can be made arbitrarily close to p by choosing N 
to grow polynomially with k. Now let us apply this to 
our model. Note that, due to the symmetry of the prob- 
lem, the number of electrons that moved to the left of the 
starting position and to the right of the starting position 
will be the same. We will therefore assume all electrons 
moving in the right direction, what can be corrected by 
an irrelevant factor of two for all length in the following 
discussion. We have a program of length l p and a quan- 
tum computer of length l q . Instead of of searching for 
l p 'electrons' to the left of the 'quantum computer' we 
can search for l p + l q 'electrons' leaving the starting area. 
The extra l q 'electrons' will be just part of the irrelevant 
code; that guarantees that the real program completely 
passed the quantum computer. Then we choose N to be 
e.g. (l p + l q ) 2 . The above calculation then tells us that 
after a time of 5000(?p + l q ) 2 the computation is success- 
ful with a probability higher than p — 0.7 (if we do not 
assume the worst case it is quite likely that we approach 
1). Since l p and l q grows only polynomial the same holds 
for the evolution time. 

Conclusion: We have shown that calculating the 
ground state energy of a translationally invariant Hamil- 
tonians in ID is as hard as solving QMA-problems. We 
have also introduced a programable quantum compu- 
tation scheme using one fixed translationally invariant 
nearest-neighbor Hamiltonian. The program is encoded 
in the initial state. The computation itself requires 
only enough patience but no active further control. We 
acknowledges support from EU projects (SCALA and 
CONQUEST), DFG-Forschungsgruppe 635, and the ex- 
cellence clusters MAP and NIM. K.V. thanks Frank Ver- 
straete for helpful discussions. 

Appendix: For equation {2} we get 



JV 



! " M 

n = NMsT, E ' 

x,y=i q,q'=i 



(q-q')(x-y)^ + At 



(3) 



where A = e(q) — e(q'). We assume that N and M are 
large and M ^> N. If we sum over the cases q — q' we get 
N/M. Since this term vanishes in the limit of M N, 
we can ignore those cases. Let us now approximate the 



sums over x and y by integrals 
dXdY cos 
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( q -q')(X-Y)N—+M 



(4) 



These integrations leads to 
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cos(At)[sin(^(g- g »] 



Nir 2 ^ 

q,q' = l 



(q - q'Y 



(5) 



Now let us define 5 = \q — q'\. Note that all terms 
only depend on the absolute value q — q' and that we 
already neglected the terms with q — q' 
get pi = jF^Y,sLi9(8)f(5). with g(5) = 



= 0. 

[sin( 1 



We 



and f(S) 



2E£i'«»(e(9) 
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e(q + S))t. Note that 
for t = we get that f(8) = 2(M — q) and in this 
case the overall sum converges to 1. One observation 
is that the main contribution to the sum comes from 
the parts where 8 is smaller than 2M/N. Furthermore 
for t > we get f(S) < 2(M - q). To get a bound 
we now want to calculate f(S) in the relevant region 
5 < 2M/N. Since 5jj is small in this region, we can ap- 
proximate e(q) — e(q + 5) by the derivative e' (q)jj lead- 
ing to f(8) = 2£^ 9 cos[sm(^)2M27r/A/]. This can 
be solved by converting the sum into an integral leading 
to f(S) = 2MJBessel(0,4^7r/M));8]. To get an upper 
bound we split up the sum in three terms 

eAI/N 2M/N M 

Pi = E 9{S)J{S) + E a(S)f(S) + J2 9(6) f (5) 

5=1 5=eM/N S=2M/N 



Assuming f(S) = 2M for the first term we get for small 
e a bound of 2e, The second term we get as bound 
uY&x[f(q)/2M,eM/N < q < 2M/N]. The last term 
we can numerical calculate for large M and get it to be 
smaller then 0.05. If we set in e = 0.001, t = 5000^ we 
get that px < 0.3. 
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